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Abstract 



■^j- ■ Recently proposed quasi-classical 9-dressing method provides a 

\ systematic approach to study the weakly dispersive limit of integrable 

^ ; systems. We apply the quasi-classical 3-dressing method to describe 

■ dispersive corrections of any order. We show how to calculate the 5 

CO problems at any order for a rather general class of integrable systems, 

presenting explicit results for the KP hierarchy case. We demonstrate 
• the stability of the method at each order. We construct an infinite set 

of commuting flows at first order which allow a description analogous 
to the zero order (purely dispersionless) case, highlighting a Whitham 
type structure. Obstacles for the construction of the higher order 
. dispersive corrections are also discussed, 

i PACS numbers: 02.301k, 02.30Jr 

a. 

1 Introduction 

During the last decade a great interest has been focused on the dispersion- 
less limit of integrable nonlinear dispersive systems. They arise in various 
contexts of physics and mathematical physics: topological field theory and 
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strings, matrix models, interface dynamics, nonlinear optics, conformal map- 
pings theory PP- [12]. Several methods and approaches have been used in 
order to study the features of this type of systems, like the quasi-classical 
Lax pair representation with its close relationship with the Whitham uni- 
versal hierarchy [3]- jjj; the hodograph transformations to calculate particular 
and interesting classes of exact solutions (rarefaction and shock waves) 0; 
the quasi-classical version of the inverse scattering method allowing to ana- 
lyze various 1+1- dimensional systems [TU]. Moreover, the recent formulation 
of the quasi-classical (^-dressing method[2U]-[2^] introduces a more general 
and systematic approach to multidimensional integrable systems with weak 
dispersion, which preserves the power of the standard method. It allows to 
build integrable systems and solutions simultaneously; that is a non triv- 
ial fact, because of the difficulty to obtain exact explicit solutions [23J. In 
addition, new interesting interrelations have been revealed, like an intrigu- 
ing connection with the quasi-conformal mapping theory, a strong similarity 
with the theory of semi-classical approximation to quantum mechanics[24 
and geometric asymptotics methods to calculate wave corrections to geo- 
metrical optics While the method has been rather widely discussed for 
purely dispersionless limit of some famous integrable hierarchies (KP, mKP, 
2DTL) [201 |2U 122 123 121] the study of the dispersive corrections is open yet. 
The latter is an interesting and challenging question in the context where 
one would like to be able to investigate the properties of the full dispersive 
system through an approximation theory. 

A main goal of the present paper is to extend the quasi-classical <9-dressing 
method to higher order dispersive corrections, starting an investigation of the 
general properties at the dispersive orders of an integrable hierarchy. We will 
show that the quasi-classical <9-dressing method works as effectively as in the 
pure dispersionless case, and the first order admits a construction parallel 
to leading order case. We present all calculations for the weakly dispersive 
KP hierarchy as illustrative example. In this paper we do not consider the 
problem of construction of bounded solutions. So, we do not care about 
secular terms which usually appear in asymptotic expansion |28j. 

We briefly present the standard <9-method in the Sec|2] and the quasi- 
classical <9-method in the SecOJ In the Sec|H we calculate the quasi-classical 
<9-problems at any order in the dispersive parameter, giving the explicit for- 
mulas at fourth order. In the Sec0we recall the method at the first order [20 , 
and prove a theorem allowing us to build an infinite set of flows reproducing 
the first corrections to the dKP equation. In the SecEJwe present the gener- 
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alization of the <9-dressing method for higher orders. The SecH contains some 
considerations about the possibility to describe the first order introducing a 
certain generalization of the Whitham hierarchy . Finally, we present some 
concluding remarks. 

2 The standard ^-dressing method. 

The standard <9-dressing method is a powerful procedure allowing to con- 
struct multidimensional integrable hierarchies of equations and their solu- 
tions. In this section we outline the main ideas of the method ic l3*U] . The 
(9-dressing method is based on the nonlocal <9-problem 

dfi A dp, x(ji, t) g(ji, t) Ro(n, fi; z, z) g~ x {z, t), (1) 

where z is the spectral parameter and z its complex conjugate, t — 
(ti,t2,tz, . . .) the vector of the times, x{ z -> z \ t) a complex-valued function 
on the complex plane C and Rq(/i, ft; z, z) is the 9-data. 

The (^-equation (0) encodes all informations about integrable hierarchies. 
It is assumed that the problem (0) is uniquely solvable. Usually, one consider 
its solutions with canonical normalization, i.e. 

+ ^ + + (2) 

z z A 

as z — > oo. 

The particular form of the function g(z,t) sets an integrable hierarchy. 
Specifically, in order to describe the KP hierarchy one has to take 

oo 

g(z,t) = e s ^ t \ S (z,t) = J2z k tk- (3) 

k=l 

Following a standard procedure 30J, it is possible to construct an infinite 
set of linear operators Mj, providing an infinite set of linear equations 

M lX = 0, (4) 

which are automatically compatible. Equations (@J) provide us with the 
corresponding integrable hierarchy. For example, the pair of operators (Lax 
pair) 
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<9y <9x 2 

d d 3 3d 3 3 _, 

M3 = — -— r U— U x 0, r U v , 

dt dx 3 2 dx 2 A x y) 



M 2 = — - — -u(x,y,t) 



with the notations 



h = x; t 2 = y; t 3 = t; 



dt 

Ub l J — oo 

provides us with the KP equation (sometimes called " dispersionfuW KP) 

d 2 u d (\d 3 u 3 du\ 3d 2 u 

+ +71TT- ( 5 ) 



dt dx dx \4dx 3 2 dx J Ady 

<9-dressing method allow us to construct a wide class of exact solutions of 
the KP equation and the whole KP hierarchy. 



3 Quasi-classical ^-dressing method. 
3.1 Dispersionless limit. 

In order to illustrate the quasi-classical d formalism, we, following Ref . |2Uj . 
introduce slow variables ti = — , d ti = ed^ and assume that g{z,T;e) = 

S {z,T) 

e i , where e is a small dispersive parameter. We are looking for solutions 
of the <9-problem in the form 

X ( z, z; — J = x(z,z; T; e) exp [ ^£i£lL I ; (6) 



e / \ e 

where x( z i z; T; e) = Y^=o fni^i T)e n is an asymptotic expansion. The 
function S has the following expansion 

a Si S 2 S 3 

S -> h^r + ^r + ...,-2-^00. (7) 

z z z z 6 
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Let us consider the <9-kernel of the following, quite general form: 



oo 

% 



Ro(fi, Jl; z,z) = - ]T(-l) fe e k ~ l - z)T k (z, z). (8) 

k=0 

The derivative 5-functions S^ k ' n ^ are defined, in the standard way, as the 
distributions such that 

-1 jf d» A d0 _ p) = (_!)*+» d ^l^\ (9) 

and := 5^°\ 
Introducing 

S(z,z;T) = S (z;T) + S(z,z;T), (10) 
one has, in the KP case, 

S = $> fe T fc + £ ^, z-oc. (11) 

k=l j=l 

A straightforward calculation for the problem (0) with the kernel (JHJ) gives 
the dispersionless (classical) limit of <9-problem 

dS „ / <9S\ 

^Ml^T > (12) 



<9z \ cfe^ 

where 

«' = E(- 1 )'(i)" r »- < 13 > 

In this limit the role of the nonlocal linear <9-problem (0) is held by the 
nonlinear <9-equation ()12j) . that is a local nonlinear equation of Hamilton- 
Jacobi type. The latter encodes all informations about the integrable sys- 
tems, like in the dispersionfull case. Nevertheless, there are substantial tech- 
nical differences in the dressing procedure. In particular, the quasi-classical 
<9-method is based crucially on the Beltrami equation's properties. We present 
them in the next subsection. 
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3.2 The Beltrami equation (BE). 

The Beltrami equation (BE) has the form 

where z G C. Under certain conditions, it has the following properties 
(see e.g.j 



1. If Cl satisfies \Q\ < k < 1, then the only solution of BE such that ^| is 
locally L p for some p > 2, and such that \I/ vanishes at some point of 
the extended plane C* is ^ = (Vekua's theorem). 

2. If \&2, ■ • ■ , are solutions of BE, a differentiable function 
f(^i, ■ ■ ■ , ^jv) with arbitrary iV is a solution of BE too. 

3. The solutions of the BE under the previous conditions give quasi- 
conformal maps with the complex dilatation Q(z. z) [3*3*j. 

It is assumed, in all our further constructions, that these properties of 
BE are satisfied. 

3.3 Zero order. 

The zero order case has been widely discussed, for different integrable hier- 
archies in various papers [201 EU 1221 HH1 1211 • Here, we outline the procedure 
for the KP hierarchy. The first crucial observation is that the symmetries of 
the (9-problem at leading order in e, for any time Tj 7 i.e. 5S = -j^rSTj, are 
given by the BE 

«(«}) (15) 

We assume that the complex dilatation W has good properties required 
for the application of the Vekua's theorem. From (fTTj) it follows that 

as j.idS! ios 2 

W 3 =z + ~zdJ) + ^dJ) + -- z "°°- (16) 

Using the BE properties, one gets the following equations in time vari- 
ables [2H1 125 
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dS 


fas 


dy 


\dx 


dS 


fdS 


~dt ~ 


\ dx 



where the notation 



is adopted and 



2 

-uo(x,y,t) = 0, (17) 

3 3 dS 
~-u —-V {x,y,t) = 0, (18) 



T 1 = x; T 2 = y; T 3 = t, (19) 



Uo = n dSl - — = -— (20) 
dx ' dx 4 dy ' 

Indeed, due to the property 2 of BE, the left hand sides of equations (fTTj) 
and (|T%|) are solutions of BE and since they vanish at z —>■ oo, then, by virtue 
of the Vekua's theorem, they vanish identically on whole complex plane. 
Equations (JT7jl and ()18j) are automatically compatible and by expansion in 
1/z-power, according to the standard <9-dressing procedure, one obtains the 
well known dispersionless KP (dKP) equation ( or Zabolotskaya-Khokhlov 
equation ) 

d 2 u 3 d f du \ 3d 2 u 

M oir- +7^- ( 21 ) 



dt dx 2 dx \ dx J 4 dy 2 

Usually the dKP equation and its dispersive corrections are obtained, by a 
slow times limit, from dispersionfull KP equation (jSJ), inserting the following 
asymptotic expansion in terms of a small dispersive parameter 

U (7) =X>nCn^- (22) 
^ ' n=0 

We stress that the formal expansion (|2*2*|) implies, in general, the appear- 
ance of secular terms. Usually one imposes additional conditions on equa- 
tions ()17j) and (j!8)l to ensure the uniform validity of this expansion [28 J. In 
this paper we do not care about boundedness of solutions (j22J) and, conse- 
quently, about secular terms. 
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Actually, as it is well known, equations ([17)1 and (|T8j) are the first ones 
of an infinite set of equations in all time variables T n , with n G N \ {0}, 
reproducing the whole KP hierarchy. Assuming 

W(z,z 1 t) = e(r-\z\)v(z,z 1 ^y (23) 

with r > 0, the function S is an analytic function outside the circle T> = 
{z G C | \z\ < r}. That is in agreement with the asymptotic behavior (jlljl . 
At last, the dKP hierarchy can be written in compact form as follows 

'z,T) = n n (p(z,T),T); n>l, (24) 



dT n 

where z G C \ V, and 



dx 

Taking into account that 



^ = ," + (i); ^oo, (27) 



then 



n n (p,T)-(^(p,T)) + = oQV z-*oo, (28) 

where 2o denotes the expansion for z obtained by inversion of equa- 
tion (J2EJ, and the symbol (•)+ means the polynomial part of the expansion. 
The left hand sides of equations ()28)1 are solutions of BE, so they vanish 
identically[22|, i.e. 



n n (p,T) = (Z$(p,T)) + . (29) 

Then, Q n can be connected to a suitable expansion series in terms of the 
p variable. Assuming 

Z (p,T)=p + jr^, (30) 

3=1 P 



8 



one has 

p\ VL 2 =p 2 + 2ai\ Viz = p 3 + 3aip + 3a 2 , (31) 
that reproduces equations (fTTj) and (fTHj). by identifications 

u = 2ai, V = 3a 2 . (32) 

The set of flows Q n represents a set of quasi-conformal maps for which the 
dispersionless integrable hierarchy describes a class of integrable deformations [21 J. 

4 Quasi-classical <5-dressing method at any 
order. 

In this section we calculate the corrections at any order for the kernel (jHJ). 

oo 

RoQi, ft z,z)= l - ]T(-l) fe e k - 1 5<»(n - z)T k (z, z). 

Let us observe that the independence of T^n on the integration variables 
\i and p,, does not reduce the generality of the kernel (jHJ). Indeed, by the 
delta function properties, it provides the same result, up to a redefinition, as 
for the kernels where I\ n depends on integration variables /i and ft. 

Substituting the expansion © into equation (JTJ), one obtains, in direct 
way, the quasiclassical <9-problem at any order. For this purpose, it is conve- 
nient to use the Faa de Bruno polynomials [3*T] defined as follows 

h n {g{x)} = d x h n -i\g(x)] + h\g(x)]h n -i\g(x)], (33) 
with h [g(x)] = 1 and hi[g(x)] = h[g(x)] = d x g(x), d x = Let us note 
that h n \g(x)] = h n [g(x)] (d x g,d^g,...,d^g\, in other words, it depends 

on the derivatives of g until n — th order. 
In our case we have 

h t [k log S(z,z)], (34) 
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for some fceN. Observing that 

i 

h l [klogS(z,z)]=Y,C l , l ,S- v , (35) 



i'=i 

we denote 



n s ,i - f so ) 



Z! 



where k is an arbitrary integer larger than I. We note that the quan- 
tity does not depend on fc. Now, we introduce the operator W p>q 



k=0 

In particular, we have 



Wp,q = — (d) ^k-p,k-q^k,n- 

h—n W 



k=0 

k 

7~Lk,k = I ~ I 



\dz J 

So, we have all ingredients to write the quasiclassical <9-problem at n — th 
order: 



(37) 

^ , (38) 

n+l 

Pn + Z^M^Wfl. (39) 
9=1 



as 


(z,z) 


- 




dz 






dz 


<9z 2 <9z 2 
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where 



It is interesting to observe that equation (|37)l is a Hamilton- Jacobi type 
equation, the first order correction is a homogeneous linear equation, while 
the next orders are linear equations, but nonhomogeneous ones. In agreement 
with the normalization condition (J2J) we have 



<P0 -> 1 + + — + — + •••> 

z Z A z 6 

, ^1,2 . ¥>1,3 . 

<Pi — + — 2- + — r + ---> 

Z Z z Z 6 

fn,l . Vn,2 <£V,3 . ^ t\t \ rnl 



The explicit problem at fourth order, i.e. n = 3 is 



D <p = 0, (41) 

A,¥>i = D m , (42) 

D Q if 2 = Dnp x + D 2 (p , (43) 

D <p 3 = Dxip 2 + D 2 ip 1 + £> 3 <p , ( 44 ) 



where 
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oz oz 2 Oz 



2 



Dl = l W "— + -W'"—- + -W" 1 — + V 4 > ( — V 

1 2 9z 2 2 ite 2 az 6 Sz 3 8 V&V 

6 <9,2 3 4 <9z 2 <9;z 2 \ 6 <9z 3 8 




24 <9z 4 12 fts 2 <9z 3 48 




24 dz 4 12 dz 2 dz 3 112 9z 3 16 V fo 2 / / 9z 2 
\2A dz* 12 dz 2 dz 3 48 
120 <9z 5 72 V ^ 3 / 48 <^ 4 4 8 V <^ 2 / <9z 3 



384 V dz2 

We stress the strong analogy between equations (|4Uj) - (J44j) and the Hamilton- 
Jacobi and higher transport equations which arise in the semi-classical approximation 
and theory of geometric asymptotics[2n]. Anyway, the main difference with 
the cited approaches is that in our case the role of the phase function is 
played by the function S that is a complex valued one rather than a real 
valued one. 



5 First order contribution. 

The BE's properties are very useful for the study of higher order corrections 
tooj201- This analysis has been initiated in the paper j20]. Some results of [201 
are presented here for completeness. Fixed an arbitrary solution S(z, z, T) 
of equation PU|). let <p and Lip be two solutions of equation (JHJ), where L 
is a suitable linear operator depending on time variables. The ratio Ltp /tp 
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satisfies the BE 



d_ f L Vo 

dz V <^ 

Then, choosing L such that 



W 



, 9 ( Lcp 



dz V 



Lip ^ 0; z — > oo 
as a result of the Vekua's theorem, one gets 

Lip Q = 0; e C. 
In particular, for the first equation of the KP hierarchy we have 



(45) 



(46) 



(47) 



L 2 <p 



A_ 2 — — - — -u(x t)) 
dy dx dx dx 2 Ul " 



d_ 

Of 



>{ dS_\ 3^ \ g ^dSd 2 S 3dS u 
dx J 2 J dx dx dx 2 2 dx 1 



where 



Ul 



t d( Po,i 



3 du 

4 cte 

9Vi 3 du i 



(48) 



<A, = 0, (49) 



(50) 



<9x <9x A dy ' 

are extracted from the condition (J46)) and % is an arbitrary solution of equa- 
tion (grp . 

Setting to zero the l/z-power expansion's coefficients in (J48|) and (J49)) . 
we find the first dispersive correction to dKP equation 



d 2 Ui 



dt dx 2 dx 2 



3d 2 , s 3d 2 Ul 

{U Ui) + 



(51) 



4 dy 2 ' 

Let us note that u\ satisfies the equation defining the symmetries of dKP 
equation (J2TJ). Indeed, considering equation (}2*T]) for uq + 5uq, one, obviously, 
gets the equation 



d 2 (du) 
dtdx 



S (Su) 



(52) 
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where 



M-) = ~M + —(■), (53) 



that coincide with ffoTj) . 

Now we will present some new results concerning the first order correc- 
tions. 

Theorem 1 Let A n and C n be the differentiable functions defined by 



A n = A n (p,T) 



d<Z8) 
dp 



1 dA 

C n = C n (p, d x p, T) = - (p, T) , 

where Z is given by B n = B n (p, T) an arbitrary differentiable 

function, ip some solution of equation and the linear operator is 
given by 

L (n) = ^ r -An-^--B n -Cn. (54) 

oT n ox 
Then L^ n '(p is also the solution of equation 

Proof. By the use of equation (}4T)j) and p = |^ (z, T) it's immediate to 
verify that 

a_K = vr aA.a_as = w ,ba. (55) 

dz dp dz dx dz 

bb. = w ,eB.e_es = w fB. 

dz dp dz dx dz 

ac. = w ,9c IL ^a_s + w ,jc^^i + w „jc^(d_a^y 

dz dp dz dx d(d x p)dx 2 d(d x p) \dz dx J 

Differentiating equation (J21J) with respect to z, we get 

d fdS\ . d (dS\ 

A n — — . (58) 



dz \dT n n dz \dx 
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Moreover, the definition of the function C„ 



c = ldAn = ldAndp + ldA L 

2 dx 2 dp dx 2 dx 



implies that 

8C n 1 dA r . 



d(d xP ) 2 dp ' ^ 

Calculating D (L^</? ), where D is given in equation (j45j). and exploit- 
ing equations (j4T) ]) .(|4*T j) and the set of equalities (}55|) - (j58|) . (j60|) . one finds 

D (i>Vo) = 0- ( 61 ) 
This completes the proof. □ 

Based on the theorem (JTJ), we choose, in particular, 

where 7] is a series defined by 

^> T ) = E^P- ( 63 ) 

3=1 P 

A choice of the coefficients in equation (|3U|) and (|63|) in such way that 
Lip —>■ for z — > oo, together with the previous arguments on the BE, gives 
us the following set of linear problems in time variables 

L<"V = 0. (64) 

These equations can be rearranged by analogy with the leading order 
case, in the form 



-£fr^(s, T) = K (q{z, T),p(z, T), d x p(z, T), T) ; 2GC\P, (65) 
where 
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Q 

Taking n = 1, 2, 3 



A„ = (c*^i I 1 dd{Z ° J ) , (66) 

\ dp 2 dx dp 

& ■■= q + rj(P, T ), 
d\ogip 



dx 



Ai = q, 

A 2 = 2pq + 2b 1 + d x p, 

3 

A 3 = 3p 2 q + 3aig + 3&ip + 3pd x p + -d x a\ + 3b 2 , 

we reproduce equations (jIH|) and P^j) and, consequently, the first order 
dispersive correction to the dKP equation by the identifications 



u = 2ai; V = 3a 2 ; 
ui = 26 x ; Fi = 36 2 . 

6 Higher order contributions. 

Unlike the first order, the higher order corrections are characterized by non- 
homogeneous equations. Here we will show how the <9-dressing procedure 
works in these cases, discussing explicitly the KP equation. 
Let us consider a set of functions <p , <p 1 , . . . ,ip n satisfying the d problem at 
(n+l)—th order and n+1 linear operators in time variables K^°\ . . . , 
such that the quantity Xlm=o ^ fn~m satisfies equation (pTTj) . Then the ra- 
tio 

Z^,m=0 Tn—m (67) 

is a solution of BE with complex dilatation W . 
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Using the same arguments as in the first order case, choosing j = 
0, . . . , n, in such a way that Ylm=o K fn-m — > for z — > oo, we get the 
linear equations 

^ (m V„- m = 0; Vz G C. (68) 

m=0 

For instance, the second order dispersive corrections to dKP equation are 
associated with the following pair 

tffV+tffVo = 0, (69) 
KfV + i^Vo = 0, (70) 

where K[ ^ = L\ and K 2 ^ = L 2 are given by (jlHj) and (jUJ), 



^ - -a?-* (71) 

' dx dx 2 \ dx 2 2 1 J dx dx 3 2 dx U2 

3dUl -V 2 , (72) 



4 <9x 



and 



0^1,1. 9V2 = 3du2 (73) 
dx ' dx A dy ' 

are deducted similarly to (J3UJ). 

Equations ()69)1 and (J70|) are compatible if and only if u 2 satisfies the 
second order dispersive correction to dKP equation 

d 2 u 2 _3g 2 3 d 2 u 2 3 d 2 uj l d 4 u 

W^~2^^ U ° U ^ + 4^ + 4~dl^ + 4~dx T ' ^ > 

The third order correction corresponds to the compatibility condition for 
the pair of the following linear problems 



K}V + ^i (1 Vi + *fVo = 0, (75) 
K [ 2 ] V 2 + AfVi + Ki\ = 0, (76) 
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where 



(2) 



(2) 



K. 



with the definitions 



-u 3 , 
d 3 



3 dS 3 d 3du, 



dx 3 2 Us dx 2 U<2 dx 



4 dx 



-V 3 



u 3 = -2 



d<p 2 ,i . dV 3 _ 3 du 3 
dx 



The equation for u 3 is of the form 

d 2 u 3 3d 2 , ,3 d 2 
{u u 3 ) + 



dx 4 dy 



. , 3d 2 u 3 ld 4 Ul 
{U1U2) + 4 



(77) 
(78) 

(79) 
(80) 



dt dx 2 <9x 2 2 <9x 2 4 <9y 2 4 dx 4 

A simple observation from equations (J69|) and (J7U|) allows us to rewrite 
them in a compact form, suggesting the generalization to any order. In order 
to realize that, we introduce the following operator 



A k [S;9] 
Noting that for k > 



e -es a k e es 

dx k 



:A k [S; 



Qk 

dx k 



S 



if 
if 



S 



k > 
k < 



r brackets 



(81) 



(82) 



the transport equations can be written as follows 



dVn E hw^ s > °K-^+i + ~ E E ^!fA fc - 2 [S; 0]^, + 



K r=0 j=0 r=0 

i=o v y 
for fc = 1, 2, 3. 

It would be nice to get a similar formula for higher times too. 



(83) 
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At last, as we noted above, a solution U\ of the first dispersive equation is 
defined up to a symmetry of the dKP equation. Using this freedom one can 
fix a "gauge" putting u\ = 0. Equations for higher corrections imply that 
the gauge u n = for each odd n is an admissible one. In such a gauge the 
even order corrective equations take the form 

d 2 iu n ) 

gtg ™ = S (u n ) + f n , n even. (84) 

That is a nonhomogeneous equation with the corresponding homogeneous 
one given by the symmetry equation. By virtue of the Fredholm's theorem, 
the nonhomogeneous term must be orthogonal to the solutions S*u of adjoint 
homogeneous equation, i.e. 

J 5*u(P)f n (P)dP = 0, P=(x,y,t,...). (85) 

One can check that for regular solutions of the KP equation, the condi- 
tion (|83j) is satisfied. Quite different situation takes place for the singular 
sector of the dKP equation. For instance, for breaking waves solutions for 
which 5u = — > oo, the condition (|83|) breaks. So, there are obstacles to 
construct the higher order quasi-classical corrections for dKP equation origi- 
nated from its own singular sector. Such an obstacle is analogous to a typical 
one appearing in the construction of global asymptotics in the semi-classical 
approximation to the quantum mechanics [24 j and in the study of the wave 
corrections to geometrical optics j2S], because of existence of caustics. In the 
dispersionless KdV case this kind of problem induces a stratification of the 
affine space of times providing a classification of the singularities |34j. The 
problem of obstacles for the construction of the higher order corrections to 
dKP equation will be considered elsewhere. 



7 A Whitham type structure. 

It is well known that the dispersionless limit of some integrable hierarchies 
admits a symplectic structure|lj-[Zj. In particular, one can see that starting 
with the function S(z,T) outside the circle V it is possible to introduce a 
2-form closed ujq in terms of the p-variable such that 

cu = dVL n (p, T) A dT n = dZ (p, T) A dM (p, T) , (86) 
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where the sum over index n is assumed. Here the last member contains 
the pair of Darboux coordinates Z G and M.q (p, T) = A4 (Z (p,T),T) = 
Aio (z, T) and the function 

M (z,T) := —(z,T) (87) 

is usually called Orlov's function 
From (IHoT). it follows that 



1 = \Z ,M \ , (88) 

J p,x 



R,2o} w , (89) 



tt n ,Mo\ • (90) 



dZ 

dT n 
dM 
dT r 

where the Poisson bracket is defined as follows 

iff a\ ( aw & f dg 9gd f fai\ 

{f( a ,0),g{ a ,0) }aJl =-----. (91) 

Equation (Jggj) is usually referred to as string equation^ and equation (|83j) 
is the dispersionless Lax pair. Moreover, equation (jSfij) implies the zero cur- 
vature condition 

lj A u = 0. (92) 

That is another way to write the compatibility condition between equa- 
tions (|89j) or (|9U|) ; indeed, expanding the total differentials, one gets from 
equation (JHU) the set of equations 

^-^ + &n,V m } P)X = 0, (93) 

that are the universal Whitham hierarchy equations)!]. In particular, for 
n = 2 and m = 3 one obtains equation (J5TJ). 

The construction of the infinite set of flows in equation (JH3J), allows us 
to describe the first correction in a way analogous to the zero curvature 
condition (|9~2*)) . Let us start with the total differential of log <p 

dlog<p {z, T) = A n {q{z, T),p(z, T), d x p{z, T), T) dT n + Mx{z, T)dz. (94) 
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The function 



Ml(z , T) = ^ {ZtT) (95 ) 

is an analogous of the Orlov's function. Differentiating equation (|94|) . one 
can introduces the 2-form 

U! := dA n {q{z, T),p(z, T), d x p(z, T), T) A dT n = dz A dMt(z, T). (96) 

The equation (j8T?j) allows us to identify a Whitham type structure defined 
by the equation 

w 1 Aa;i = 0, (97) 
which can be given explicitly 



^A + (A A] +IAAI — + 1 A A i ^t- 
= G n (A m )-G m (A n ), (98) 
where G n is defined as 

Let us note that for n = 2 and m = 3, equation (|98|) gives, of course 
equation (j5T|) . 



8 Concluding remarks and perspectives. 

In this paper we demonstrated how to generalize the quasi-classical <9-dressing 
method at any dispersive order. In particular, the formula (J65|) suggests that 
a regular structure survives at higher orders, allowing to define a potential en- 
coding all informations about the hierarchy at the first dispersive order, just 
like in the purely dispersionless case. Moreover the paper provides us with 
several intriguing observations which could be subjects of future study. In 
particular, it would be useful to analyze the connection between the standard 
theory of the asymptotics approximation and the quasi-classical (^-dressing 
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approach. This could set a relationship between the source of the obstacles in 
the construction of the corrections, discussed in the section ©, and the caus- 
tics problem. As far as concerning possible applications it would be of interest 
to consider the KdV limit of the KP hierarchy to establish a connection with 
the Dubrovin- Zhang theory [3]. Another intriguing matter of study and appli- 
cation is associated with a possible physical interpretation of the solutions of 
dispersionless systems as describing integrable dynamics of interfaces ^7]. In 
this context, one should investigate a quasiconformal dynamics described by 
the quasiclassical 5 problem and the corrections (|37H39|) . One could analyze 
the singular behavior of interfaces' evolution, by introducing, for instance, 
a small "dispersive" parameter. For such a purpose, a deeper study of the 
singular sector order by order will be also necessary. 
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